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ABSTRACT 

There is growing observational evidence that the variability of red giants could be 
caused by self-excitation of global modes of oscillation. The most recent evidence of 
such oscillations was reported for a UMa by Buzasi et al. (2000) who analysed space 
photometric data from the WIRE satellite. 

Little is understood about the oscillation properties in red giants. In this paper 
we address the question as to whether excited radial and nonradial modes can explain 
the observed variability in red giants. In particular, we present the results of numerical 
computations of oscillation properties of a model of a UMa and of several models of a 
2M star in the red-giant phase. 

The red giant stars that we have studied have two cavities that can support 
oscillations: an inner core that supports gravity (g) waves and a surrounding shell that 
supports acoustic (p) waves. Most of the modes in the g-mode frequency range are g 
modes confined in the core; those modes whose frequencies are close to a corresponding 
characteristic frequency of a p mode in the outer cavity are of mixed character and 
have substantial amplitudes in the outer cavity. We have shown that such modes of low 
degree, I = 1 and 2, together with the radial (p) modes, can be unstable. The linear 
growth rates of these nonradial modes are similar to those of corresponding radial 
modes. In the model of a UMa and in the 2M Q models in the lower regions of the 
giant branch, high amplitudes in the p-mode cavity arise only for modes with £ = 2. 

We have been unable to explain the observed oscillation properties of a UMa, 
either in terms of mode instability or in terms of stochastic excitation by turbulent 
convection. Modes with the lowest frequencies, which exhibit the largest amplitudes 
and may correspond to the first three radial modes, are computed to be unstable if all 
effects of convection are neglected in the stability analyses. However, if the Lagrangian 
perturbation of the turbulent fluxes (heat and momentum) are taken into account in 
the pulsation calculation, only modes with higher frequencies are found to be unstable. 
The observed frequency dependence of amplitudes reported by Buzasi et al. (2000) does 
not agree with what one expects from stochastic excitation. This mechanism predicts 
an amplitude of the fundamental mode about two orders of magnitude smaller than 
the amplitudes of modes with orders n > 5, which is in stark disagreement with the 
observations. 
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1 INTRODUCTION 

Buzasi et al. (2000) have reported the discovery of oscilla- 
tions in photometric data from a UMa observed with the 
star camera on the WIRE satellite. They interpret the ob- 
served oscillations as radial modes, and cautiously suggest 
that the modes may be excited by the mechanism similar to 
that responsible for solar oscillations. The star, however, is 
very different from the Sun. Its spectral type is K0 III. Mod- 



els of the star's internal structure, and its pulsation proper- 
ties, suggest that the star is a red giant. Thus, even if the 
oscillations are stochastically excited by turbulence in the 
outer convective zone, as they are in the Sun, some impor- 
tant differences between the oscillation properties of a UMa 
and the Sun should be expected. 

Variability is a common feature of red giants. There 
are strong observationally based arguments that, at least in 
part, this variability is due to global pulsations. Edmonds 
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& Gilliland (1995) proposed radial or nonradial pulsations 
as an explanation for the variability they observed in K gi- 
ants in the globular cluster 47 Tuc. They found frequen- 
cies between 3 and 6 /iHz with amplitudes between 5 and 
15 mmag. Cook et al. (1997) analysed photometric data from 
red stars in the LMC collected from the microlensing project 
MACHO, and reported that the period-luminosity relation 
comprises several ridges in the period range of 10 - 200-days, 
which may be interpreted as arising from radial modes of 
different n. Hatzes and Cochran (1998) reported that there 
is strong evidence from radial-velocity data for oscillations 
in K giants. Frequencies similar to those found in a UMa 
have been found in a number of other objects. In particular, 
Hazes and Cochran quote 11 frequencies of Arcturus (K2 III) 
in the range 1.4-6.8/iHz. Evidence for short-period multi- 
mode pulsations in a number of M stars has been presented 
recently by Koen and Laney (2000) . 

Thus, by being a multimode pulsator a UMa appears 
not to be unique amongst red giants. But with its high fre- 
quencies and low amplitudes it does represent an extreme 
case so far, although it is not wholly out of line with the 
others. This star is the hottest and the least luminous ob- 
ject amongst variable red giants. It provides, so far, the best 
example of possible red-giant oscillations, but its spectrum 
is not as clean as that of the Sun, or of those of other main- 
sequence stars, or of white dwarfs. Nevertheless we consider 
the evidence to be strong enough to justify new investiga- 
tions in the theory of red-giant oscillations. So far, only the 
modelling of radial pulsations in Miras (see Xiong et al. 1998 
and references therein) and Arcturus (Balmforth et al. 1991) 
has attracted the theorists' attention. Nonradial oscillations 
in red giants have been ignored almost entirely. Here we re- 
view the theoretical aspects of this problem, in Section 3, 
and provide some numerical examples of the properties of 
the oscillations of a model of a UMa and of models of a 
2Mq star on the red-giant branch. Some data concerning 
these models are presented in Section 2. 

The most intriguing issue posed by the discovery of os- 
cillations in red giants is the identification of the mecha- 
nism by which they are driven. Possibilities to consider are 
(a) stochastic excitation of linearly stable modes by con- 
vection and (b) self-excitation of linearly unstable modes. 
We shall speak of oscillations of case (a) as being solar-like, 
and case (b), Mira-like. Our understanding of the excita- 
tion mechanism in the Sun and in Miras is not satisfactory, 
but the separate association of these stars with each of the 
two distinct possible excitation mechanisms is now generally 
accepted. In Section 4 we present results of calculations of 
radial-mode stability and of the amplitudes in the case of 
stochastic excitation. 



2 SELECTED MODELS 

There are stringent constraints on the parameters for defin- 
ing models of a UMa. The star is bright, and is in a visual 
binary system. Accurate spectroscopic data, parallax, and a 
radius determination by means of interferometry are avail- 
able. After considering all the observational data, Guenther 
et al. (2000) suggested the following values for the star's 
global parameters: M=4-5M©, log(L/L Q ) = 2.5 ± 0.05, 
r off = (4660 ± 100) K. Guenther etal. (2000) constructed 
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Figure 1. H-R diagram showing evolutionary tracks of models 
with masses of 2 and 4 Mq. Models selected for pulsation anal- 
yses are indicated by diamonds. The box around the symbol for 
the 4 star indicates the uncertainty in locating a UMa in the 
diagram (see Table 1). 

evolutionary models with masses in this range and with an 
initial chemical composition X — 0.727 and Z = 0.0124, 
which is consistent with the spectroscopic value of [Fe/H] 
and the Galactic helium enrichment. They found that only 
models with M < 4.5Mq satisfy the observational con- 
straints. 

We have adopted the same initial chemical composition 
in our model calculations. Furthermore, we have adopted the 
same opacity and equation of state. For the model of a UMa 
we have considered only M = 4Mq, and we have adjusted 
the mixing-length parameter, a, to be consistent with the 
values of log(L/L Q ) and T e s proposed by Guenther etal. 
(2000). 

The star a UMa is a high-mass red giant with a non- 
degenerate core. Such stars are very rare. As seen in Fig. 1, 
the red-giant branch for M — 4 Mq is very short; the star 
spends only 0.5 My on it, which is more than two orders of 
magnitude shorter than the time spent by a star with a mass 
of M = 2Mq, in which helium ignites in a degenerate core. 
We have chosen the model sequence with M = 2 M to il- 
lustrate nonradial mode properties in red giants over a wide 
range of luminosity. The most important parameter deter- 
mining nonradial mode properties is the ratio of the mean 
density of the core to the mean density of the whole star. 
In the sequence we have chosen, this parameter increases by 
nearly four orders of magnitudes between the bottom and 
the top of the giant branch. We have considered four models 
for the 2 Mq sequence, calculated with the same values of 
X, Z and a as those for the model of a UMa. The locations 
of the selected models on the evolutionary tracks are indi- 
cated in Fig. 1; the parameters characterizing these models 
are listed in Table 1. 



3 NONRADIAL MODES OF RED GIANTS 

Guenther et al. (2000) considered only radial modes as po- 
tential candidates for explaining the peaks in the a UMa fre- 
quency spectrum determined by Buzasi et al. (2000) . They 
noticed that the frequencies of these peaks are much lower 
than the buoyancy frequency deep in the star, and presumed 
that a nonradial interpretation would appear to imply that 
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Table 1. Some parameters of the models used in this work; M c is the core mass including the hydrogen burning shell; the subscript be 
denotes bottom of the convective envelope. 



Model 


M/Mq 


age(Gy) 


logTeff 


log L/L Q 


R/Rq 


logT c 


logpc 


M c 


rbc/R 


M hc /M 


M a 


4 


0.143 


3.6993 


2.50 


27.21 


7.924 


4.187 


0.1083 


0.3259 


0.4300 


M21 


2 


0.904 


3.6915 


1.50 


7.77 


7.753 


4.766 


0.1145 


0.1481 


0.1869 


M22 


2 


0.934 


3.6596 


2.00 


15.99 


7.776 


5.220 


0.1397 


0.0509 


0.1527 


M23 


2 


0.963 


3.6240 


2.50 


33.50 


7.780 


5.555 


0.1705 


0.0264 


0.1755 


M24 


2 


0.972 


3.5848 


3.00 


71.32 


7.835 


5.752 


0.2017 


0.0124 


0.2042 



the modes are g modes of high radial order. According to 
their estimate, the separation between the cyclic frequen- 
cies of consecutive g modes of like degree is of the order 
of 0.1 fiRz, and consequently, they argued, the spectrum 
could not be resolved into individual modes. They did not, 
however, explain why radial modes should stand above this 
quasi-continuum, which would appear to be necessary for 
explaining Buzasi's observations. Moreover, they failed to 
point out that a g mode that resonates at a corresponding 
(i.e., same value of £) characteristic p-mode frequency of the 
outer acoustic cavity can have a particularly large amplitude 
at the surface. 



3.1 General properties 

The basic properties of nonradial oscillations in highly 
evolved stars were determined in the 1970s (Dzicmbowski, 
1971, 1977; Osaki, 1977). However, the objects of interest in 
those early works were stars in the Cepheid instability strip. 
To the best of our knowledge there is only one paper devoted 
to the theory of nonradial oscillations in red giants. It is a 
short note by Keeley (1980), in which a crude estimate of 
mode trapping was made. The nonadiabatic effects, which 
are very important in this context, were ignored. 

The differences in the nonradial mode properties be- 
tween red and yellow giants are a consequence of the differ- 
ent depths of the convection zones. The formalism for cal- 
culating linear modes in these two types of star is the same. 
Here we provide only an outline of the formalism developed 
by Dziembowski (1977), which was recently recalled in some 
detail by Van Hoolst et al. (1998, hereafter VDK). We intend 
to apply this formalism to low-degree modes (I = 1,2), in 
a cyclic- frequency range starting somewhat below the fun- 
damental radial-mode frequency and extending up to the 
acoustic cut-off frequency, ^ ac , in the photosphere (r = R). 
In our model of a UMa, v ac — 27pMz. Buzasi et al. (2000) re- 
ported peaks in the power spectrum located above our value 
of ^ ac , which evidently cannot easily be interpreted in terms 
of strongly trapped acoustic modes. 

The starting point of our discussion is an asymptotic 
solution, for large order n, of the nonadiabatic wave equa- 
tion, which is valid in the radiative interior. In this approx- 
imation, any perturbed scalar parameter may be expressed 
with respect to spherical polar coordinated (r, 9, 0) in the 
following form: 



q(r,6,4>,t)=A(r) [e 1 *^ + e~ 1 *^] Y™(9, <f>)e 



iuj c t 



(i) 



in which t is time. The amplitude, A, is a slowly varying 
function of r. The rapid variations are described by the 



phase $, which for stars with radiative cores may be written 
in the form 



<E>(r) 



Jo 



k dr — 



'+1W. 



(2) 



The general expression for the radial component, k, of the 
wave number of high-order modes in the gravity-wave cavity 
may be found in VDK. The quantity lo c = ui — i~f (with lu > 
0) is the complex eigenfrequency. We focus our attention on 
predominantly oscillatory modes, i.e. modes with a growth 
rate 7 satisfying the condition 7I <g; lu. We also specify q 
to be the relative Lagrangian perturbation to the pressure, 
Sp/p. 

If the radiative energy losses are regarded as being 
small, there is a simple expression for the radial wave num- 
ber far from the edges of the cavity: 



lu r lu 

where 



N 




(3) 



(4) 



is the buoyancy (Brunt- Vaisala) frequency, g is the 
local gravitational acceleration, p is density, Ti = 
(d In p/d\n p) a( j is the first adiabatic exponent, and 



l(l+l) fl L r V ad 



(5) 



in which L T (r) is the total rate at which radiant energy 
crosses a sphere of radius r, V = din T/dln p and V a d = 
(<9ln T/dln p) a a- The quantity T> is a measure of the radia- 
tive energy loss, which is the only nonadiabatic effect we 
consider in this cavity. Here, V ac i > V is always satisfied; 
hence V > 0. To the same approximation, the amplitude is 
given by 



A(r) oc 



V 



(6) 



r 3 \fkp~ 

with V — grp/p. The approximate expressions for the wave 
number (equation 3) and for the amplitude (equation 6) 
are valid only if V < 1. For models M23 and M24 (see 
Table 1) the computations suggest that T> ^> 1 in certain 
layers inside the star, at least for some modes considered 
in the calculations. However, for these models we use this 
approximation only in the outer layers of the asymptotic re- 
gion, where the approximation is satisfied. The maximum 
value of V in our model of a UMa for the lowest-frequency 
quadrupole (£ = 2) mode is 0.2. As in all red-giant models, 
that maximum occurs within the shell source. Even if V is 
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small, consequences of radiative losses may still be very im- 
portant for the wave properties, because if N/cu is large we 
may still have $i = Ss(<E>) >• 1. Equation (1) describes a su- 
perposition of an outward (first term on the rhs) and inward 
propagating gravity wave, where the direction of propaga- 
tion is the direction of the group velocity. 

Let us concentrate now on the oscillations in the outer, 
acoustic cavity. Moreover, let us assume that the approxi- 
mation for q given by equation (1) is valid in an interval 
[0, T{\ of the g-mode cavity. In this interval we can neglect 
the derivative of A with respect to r in the calculation of 
the derivative of q. Thus, at r = n we have approximately 



ldq 
q dr 



ik 



cxp(i$) — exp(— i$) 
cxp(i$) + cxp(— i$) 



(7) 



which provides a boundary condition for the numerical solu- 
tion of the equation for the linear nonadiabatic oscillations 
in the interval [rf,i?], which contains the very outer layers 
of the g-mode cavity in which the asymptotics breaks down, 
the entire surrounding acoustic cavity and the evanescent 
zone between. The lhs of equation (7) depends on I and on 
the value of uj c . The dependence on lo c is relatively weak 
in comparison with the explicit dependence of the rhs of 
equation (7). 

Assuming <3>; 3> 1, expression (7) simplifies to 



8q 
dr 



-\kq, 



(8) 



which is valid for the case when the wave is effectively dis- 
sipated on its way towards the centre. The energy loss may 
be overcompensated by the driving operating in the outer 
layers if the wave amplitude A(rf) is small, i.e., if the mode 
trapping in the acoustic propagation zone is severe. Indeed, 
nonradial modes with growth rates 7 similar to those of ra- 
dial modes were found in models of Cepheids (Osaki, 1977) 
and of RR Lyrae (Dziembowski, 1977). In these two indepen- 
dent papers, equation (8) was used for the inner boundary 
condition. Such modes were named in VDK as S(trongly) 
T(rapped) U(nstable). We must emphasize, however, that 
even when the amplitude in the outer acoustic zone is rel- 
atively large, according to VDK the oscillations typically 
have 80% of their energy in the inner, g-mode cavity. In 
the model of the RR Lyrae star considered by VDK, STU 
modes were found only with £ > 4. We shall see that in red 
giants STU modes may exist also with I as low as unity. 
The frequency separation between consecutive low-degree 
STU modes is similar to the separation between consecutive 
radial modes. 

The STU modes are true eigensolutions of the nonadi- 
abatic oscillation equations for the whole star. Indeed, the 
boundary conditions (7) and (8) are equivalent for unstable 
modes because if 7 > we have Q(k) > throughout the 
interval [0, n], and $i does not change sign. If V <C 1, this 
conclusion follows immediately from equation (3) although, 
in fact, it is true also for any T> > (see e.g. Dziembowski, 
1977). For stable modes, the situation is more involved. If a 
solution with 7 < is found subject to the inner boundary 
condition (8), then the solution must always be checked to 
determine whether it satisfies the inequality $i ^> 1. Actu- 
ally, this inequality is rarely satisfied. Let us note that with 
the usage of equation (7) we assume maximum energy losses. 
When we use equation (7) with a properly calculated phase, 



instead of equation (8), we may find unstable modes. How- 
ever, for such modes the growth rates are typically much 
smaller than those of radial modes with similar frequencies. 

A dense spectrum of weakly unstable modes with £ = 1 
and 2 was found for the RR Lyrae model considered by 
VDK. We shall discuss in the next Section the problem of 
mode stability, and we shall see that it is actually far from 
being solved. Fortunately, whatever the mechanism respon- 
sible for the excitation of the modes, it should operate in 
the layers where the radial eigenfunctions do not depend on 
I. We shall take advantage of this property in our discussion 
of the relative chances of nonradial or radial modes being 
excited. 

It seems to be not unreasonable to assume that if in a 
certain frequency range unstable modes of various degrees 
exist, their chances of being excited are related to the growth 
rates 7. The growth rate may be expressed in terms of the 
work integral, W , and the mode inertia, J, through the well- 
known relation (see e.g. Unno et al. 1989) 



7 : 



W 
2luI' 



(9) 



The generic expression for the work integral is 



W 



J d 3 xp[ 



-TV ad 3(<f <5s)] 



6 ^ 5 -^) ,(10) 

P Pt J 



where s is the entropy per unit mass and p t (the so-called 
turbulent pressure) is the rr-component of the Reynolds 
stress tensor Tij = pmuj (u is the turbulent velocity field 
and the overbar denotes an ensemble average). The asterisk 
denotes complex conjugate. In this expression we neglect the 
contribution from the anisotropy of the Reynolds stress ten- 
sor, which is small compared with the isotropic component. 
But in this section we neglect convection dynamics in the 
model computations: in the evaluation of the work integral 
the second term of the rhs of equation (10) is neglected. The 
mode inertia is defined as 



/ 



d 3 xp||j 3 



(11) 



with £ representing the displacement eigenfunction. The in- 
tegrals are over the entire volume of the star. The inertia / 
enters also into the expression for the amplitudes of stochas- 
tically excited modes (see equation 21). Another quantity in 
the expression for the amplitudes is the energy supply rate 
Pq injected into the modes by the turbulent convection, and 
which we assume to be generated predominantly by the fluc- 
tuating Reynolds stresses (see next Section). 

There are important differences between radial and non- 
radial oscillation properties below the acoustic propagation 
zone of the nonradial modes. These differences are reflected 
in the values of /. If / is large, the largest contribution to 
the work integral may arise in the gravity-mode propagation 
zone, where the asymptotic approximation is applicable. In 
the gravity-wave propagation zone we have adopted for the 
Lagrangian specific entropy perturbation the expression 



8s = 2ic v 



dlnT 
9 hip 



dlnp 1 
dlnp ~ fT 



Vq, 



(12) 



in which c p is the specific heat at constant pressure. It fol- 
lows, for example, from equation (19) of VDK in the weakly 
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nonadiabatic limit, that is, when \V\ <gC 1. Substituting 
equations (1) and (12) into equation (10) we obtain 



W g = -2CW Z 



hVd$ r , 



(13) 



where 3> r = 5ft(<3>) and C{ is a real positive constant which is 
obtained from the eigenfunctions calculated numerically for 
r > n, and 



h = exp(2$i) + exp(-2$i) . 



(14) 



The oscillatory term, proportional to cos(2$ r ), in the inte- 
grand of W g has been ignored, which is consistent with the 
asymptotic approximation. 

An expression for £ in terms of q can be calculated in 
the adiabatic approximation. The result is 



? V 



/e r + 



t(t + 1) 



V h 



(15) 



where e r is a unit vector in the radial direction and Vh is 
the horizontal component of the gradient operator. From 
equation (3) we conclude that, to a first approximation, the 
contribution of the radial displacement to I can be neglected 
if N 3> u>. Thus, from the asymptotic interior we obtain the 
following contribution to the modal inertia / : 



J o 



(16) 



We normalize the relative rms radial component of the 
mode displacement to unity at the stellar surface. The coeffi- 
cient C{ , which is a function of w, exhibits minima separated 
by nearly the same interval in frequency as the minima in the 
radial modes. This is a manifestation of the trapping prop- 
erties of the acoustic cavity, which are purely dynamical and 
result from a resonance between the two cavities. The spec- 
trum of g modes in the inner cavity is so dense that there is 
always a g-mode-like oscillation whose frequency resonates 
with a p mode in the outer cavity, such that the amplitudes 
in both cavities are similar. All other g modes are confined to 
the central g-mode cavity, and have very low amplitudes in 
the outer layers of the star. Mode trapping is influenced also 
by the behaviour of the factor h (see equation 14), which is 
determined by nonadiabatic effects. For STU modes we have 
h « exp(2<E>i,i) 3> 1, which is a sharply increasing function 
of r. Thus, I g is negligible and W g may be evaluated as the 
rate of wave losses: 



W R 



''I 



%{£q) sin(9d(9d0 



-Cuo hf . 



(17) 



For all other modes we have to use equations (13) and (16) to 
evaluate the contributions W g and I g . Equation (3) implies 
that for stable modes h(r) has a maximum in the layer in the 
star in which T> — —7 /to. Thus, 7 g may be a significant, and 
is often the dominant contribution to /. Let us note that the 
values of 7 g and W g depend on 7, and consequently on the 
nonadiabatic processes operating in these layers. This means 
that uncertainties in the computation of the nonadiabatic 
effects are to some degree reflected in the values of I g and 
W g . Damping in the outer layers reduces the effect of mode 
trapping. 




10 14 1 i 

v GuHz) 

Figure 2. Modal inertia in units of 3MR 2 , plotted against fre- 
quency. The eigenfunctions arc normalized such that at the sur- 
face £ r = RY l m (8, (p) cxp(iw c t). Individual nonradial modes are 
not resolved. The symbols are displayed only for those modes 
that are locally most trapped. D p = means that all nonadia- 
batic effects in the outer layers are ignored. 



3.2 Application to a UMa 

In our code for computing nonradial nonadiabatic oscilla- 
tions (Dziembowski, 1977) we set the Lagrangian perturba- 
tion of the turbulent fluxes (heat and momentum) to zero, 
and we ignore the turbulent pressure in the equilibrium 
model. With this treatment, all radial modes are found to 
be unstable. 

In Fig. 2 we show the behaviour of the normalized mode 
inertia J n = I/3MR 2 as a function of the cyclic frequency 
v — ui/2ty for our model of a UMa. The choice of normal- 
ization is not important here, except that all modes are as- 
sumed to have the same surface amplitude of radial displace- 
ment. There are two sequences of model results: in the first 
sequence (upper plots) we calculated h with 7 obtained from 
our code; in the second sequence we suppressed all nonadia- 
batic effects where r > n . A comparison allows us to assess 
some of the consequences of the uncertainties of the physics 
in the convective zone. 

Symbols are used to denote the nonradial modes that 
are most trapped in the acoustic cavity. For the I — 2 se- 
quence the minima in 7 n almost coincide with the radial- 
mode frequencies, while the minima for I = 1 are located 
roughly half-way between the i = and 1 = 2 minima. The 
positions of these minima resemble the positions of modes in 
the whole-disc spectra of solar oscillations. There are many 
more modes with I = 1 and I = 2 than are depicted by the 



© 0000 RAS, MNRAS 000, 000-000 



6 W. A. Dziembowski, D. 0. Gough, G. Houdek, R. Sienkiewicz 





Figure 3. Contributions to the energy dissipation rate from the 
g-modc propagation zone in units of the stellar luminosity, L . In 
the upper panel the energy gain rate for radial modes is plotted 
with open circles. See caption of Fig. 2 for further information. 






i 
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vy/rx/G<p> 

Figure 4. Occurrence of instability in the sequence of 2 Mq gi- 
ants of nonradial mode of low degree I and low radial order n 
that are most strongly trapped in the outer acoustic cavity. The 
abscissa is the dimensionless frequency y^ir/G < p >u. Full sym- 
bols denote unstable modes, open symbols stable modes. Circles, 
triangles, and squares denote I = 0,1, 2 modes, respectively 



When I = 2 the modes that are most trapped are de- 
tached from the remaining modes, except for the one at 
v ~ 10/iHz. Except for this particular mode, all the other 
modes satisfy $i,f > 1. Thus, the unstable 1 = 2 modes are 
STU modes, and their growth rates are nearly the same as 
those of the corresponding radial modes. All the unstable 
I = 1 modes have <3?i,f <C 1, and the trapping effect is less 
severe. The inertiae of the most trapped I = 1 modes are 
always significantly larger than those of the closest radial 
mode (see Fig. 2). 



symbols. The frequency separation Av between nonradial 
modes of consecutive order n is indeed very small. It may 
be evaluated from the asymptotic formula 

^=2.4x10 4 - v (18) 

where v is expressed in /iHz. The numerical constant is spe- 
cific to the model. At I = 1 and v = 2.8piHz the value of Av 
is about 0.0013 /iHz, much less than that found by Guenther 
etal. (2000). 

There is no substantial difference between the trapping 
pattern of the two sequences, except for the differences in 
the depths of the minima, particularly those of the I = 2 
modes. Greater driving in the outer layers results in deeper 
minima. If there is net damping in the outer layers, as in 
the case of solar oscillations, the minima are shallower than 
in the adiabatic approximation. 

In Fig. 3, we plot the rate of energy dissipation, D s = 
—LoWg, in the asymptotic interior for the same two sequences 
of modes. In addition, in the upper panel, we show the total 
energy gain rate, — D p = —ujW, for radial modes. The total 
driving rate for the nonradial modes is given approximately 
by 7 ~ (D p (v) + D g ) /2ui 2 1 , because significant nonadiabatic 
effects may arise only either near the surface - and then they 
are I- independent - or in the deep interior where the g-modc 
asymptotics applies Some nonradial modes with frequencies 
larger then 12 /iHz are found to be unstable: if the radial 
modes with v > 12 /iHz are indeed unstable, then there are 
also some unstable low-degree nonradial modes. 



3.3 "Unstable" low-degree modes in 2 M© red 
giants 

In Table 2 we compare some characteristics of modes of the 
a UMa model, M a , and the modes of models of 2 M Q red 
giants that are found to be unstable with our code. For the 
models M21 and M22 we find more-or-less similar properties 
to those of the M a model. Strong trapping occurs only for 
£ > 2, and the nonradial modes are unstable for the higher 
n. Note that n is the radial order only for £ = (for £ > 
it is the number of nodes in the acoustic cavity of the radial 
component of the displacement eigenfunction) . For all cases 
the upper limit of the unstable range is determined by the 
acoustic cut-off frequency. 

In the more luminous giants (models M23 and M24) 
STU modes are found even for I = 1. In Fig. 4, we show 
how instability of the most strongly trapped modes increases 
with stellar luminosity. 

The relative frequencies of the most strongly trapped 
modes of the models considered here are different from those 
of the RR Lyrae star model considered by VDK and of 
RR Lyrae stars in general. In red giants the most strongly 
trapped £ = 1 modes are located between the £ = and 
1 = 2 pairs, whose frequencies are nearly coincident. There 
is a similarity with solar p modes, although in the case we 
have studied here the frequencies of the strongly trapped 
1=1 modes are somewhat closer to the higher-frequency 
even-degree pair. In RR Lyrae stars, on the other hand, the 
frequencies of the most strongly trapped £ = 1 modes are 
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Table 2. "Unstable" low-degree modes in red-giant models; u ac is the cyclic acoustic cut-off frequency computed for an Eddington grey 
atmosphere; IIo,i is the period of the fundamental radial mode; for £ > 0, n indicates the range of consecutive modes that are most 
trapped. 



Model 




n ,i 




n-range 






f-rangc (/iHz) 






/iHz 


days 


i = o 


1=1 


1 = 2 


e = o 


1 = 1 


t = 2 


M a 


26.6 


4.10 


1-14 


8-14 


4-12 


2.82-26.3 


16.2-27.7 


15.0-28.2 


M21 


178. 


.919 


1-18 


6-18 


6-18 


12.6-166. 


59.5-162. 


62.5-166. 


M22 


43.4 


2.65 


1-13 


5-14 


3-13 


4.36-40.3 


16.3-41.9 


11.5-40.3 


M23 


10.2 


7.74 


1-9 


2-10 


1-12 


1.50-9.30 


2.37-9.81 


1.49-10.1 


M24 


2.25 


24.1 


1-6 


1-6 


1-6 


.480-2.03 


.397-1.92 


.483-2.01 



actually closer to the radial eigenfrequencies than are the 
frequencies of the most strongly trapped 1 = 2 modes. 

These differences between RR Lyrae stars and red gi- 
ants are related to red giants having much deeper convec- 
tive zones. The differences are reflected in the different be- 
haviours of the Brunt- Vaisala frequency, N. In Fig. 5 wc 
compare TV and the Lamb frequencies L\ and L2 in the en- 
velope of an RR Lyrae model with those in the envelope of 
model M21. A deeper convective zone is associated with a 
wider evanescent zone separating the p-mode and g-modc 
propagation zones, and hence there is a possibility of more 
efficient trapping. This is why we find STU modes of low 
degree in red-giant and not in RR Lyrae models. Exception- 
ally poor trapping of the £ = 2 modes in the frequency range 
of the first two radial modes is due to the narrowness of the 
evanescent zone. 



4 EXCITATION MECHANISMS 

There is little doubt that the interaction between pulsation 
and convection plays an essential role in red-giant pulsation, 
and that the approach adopted by us to obtain the results 
reported in the previous sections is inadequate. The driv- 
ing agent that caused instability of the radial modes is the 
same as that suggested first by Ando and Osaki (1975) in 
an attempt to explain solar p-mode excitation in the Sun, 
and is artificial. It is easy to understand why: at the pho- 
tosphere, where the energy is carried mostly by radiation, 
the flux perturbation is negative in the high-temperature 
phase of the pulsation cycle. This is a result of the steep 
increase of the opacity with temperature in the outer layers. 
The fraction of the energy carried by convection increases 
rapidly inwards. Since, by assumption, the convective flux re- 
mains unperturbed, the energy is forced to be captured by 
the photospheric layers, and the putative heat engine works. 
This phenomenon is sometimes called convective blocking, 
which is confusing because what actually blocks the heat 
flux is the opacity variation. However, there is no physi- 
cal justification for the neglect of the perturbed convective 
heat flux and Reynolds stresses. Indeed, pulsational modu- 
lation of the convectively unstable stratification of the star 
is bound to modulate the convective dynamics, and domi- 
nate the driving or damping in regions where the convective 
fluxes dominate in the equilibrium state. 

Effects of convection on the stability of radial pulsa- 
tions in cool stars have been investigated since the early 
1970s (see e.g. Xiong etal. 1998; Houdek, 2000). Recent ef- 
forts have focused mainly on Mira stars and the Sun. Ac- 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



r/R 

Figure 5. Lamb frequencies L( = yJJQ. + T)c/r (for £ = 1 and 
2) and the Brunt- Vaisala frequency N in the M21 model and in a 
representative RR Lyrae star model. The latter is characterized by 
the following parameters: M = O.67M , Y = 0.243, Z = 0.001, 
Y c = 0.17 (helium abundance in the core), log(L/L0) = 1.717, 
logT e ff = 3.822. The ordinate scale is dimcnsionlcss, and corre- 
sponds to angular frequencies measured in units of ^4nG< p>, 
which corresponds to the dimensionless cyclic frequencies of 
Fig. 4. In these units the frequencies of the first two radial modes 
in the RR Lyrae model are 1.82 and 2.46. 

cording to the calculations of Xiong etal. (1998), low-order 
radial modes of Mira models are unstable, whereas those of 
orders n > 4 were always found to be damped (see also the 
work by Balmforth etal. (1991) on Arcturus). In a study 
of p-mode stability in the Sun by Balmforth (1992a), all 
modes have been found to be stable. Balmforth used in his 
calculations Gough's (1976, 1977) nonlocal, time-dependent 
mixing-length model for convection, improving on the code 
used by Baker & Gough (1979) to study RR Lyrae stars 
by incorporating the Eddington approximation to radiative 
transfer for both the equilibrium structure and the pul- 
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sations. Houdek etal. (1999) applied these calculations to 
solar-type stars, and estimated amplitudes of intrinsically 
stable stochastically excited radial oscillations in stars with 
masses between 0.9 M and 2.0 M close to the main se- 
quence. 



4.1 Linear stability of radial modes in aUMa 

Here we apply Balmforth's (1992a) treatment of pulsation 
to a model of a UMa. In particular, we include turbulent 
pressure in the equilibrium model, and the stability analy- 
sis includes the Lagrangian perturbations of the convective 
heat and momentum fluxes. We use an envelope model cal- 
culated with the surface parameters of model M a given in 
Table 1, and an atmosphere using the T-t relation of model 
C of Vernazza, Avrett & Loeser (1981). The value of the 
mixing-length parameter was adjusted such as to reproduce 
the same depth of the convective zone as was obtained from 
the evolutionary computation. The nonlocal treatment of 
convection introduces two more parameters, a and b, which 
characterize respectively the spatial coherence of the ensem- 
ble of eddies contributing to the total heat and momentum 
fluxes and the extent over which the turbulent eddies ex- 
perience an average of the local stratification. Theory sug- 
gests approximate values for these parameters, but it is ar- 
guably better to treat them as free. Roughly speaking, the 
parameters control the degree of 'nonlocality' of convection; 
low values imply highly nonlocal solutions, and in the limit 

a, b — > oo the system of equations reduces to the local formu- 
lation (except near the boundaries of the convection zone, 
where the local equations are singular). 

The energy dissipation rate D p of radial p modes was 
calculated as a continuous function of oscillation frequency 
by relaxing the inner dynamical boundary condition. The 
results shown in the upper panel of Fig. 6 were obtained 
for two sets of the nonlocal convection parameters a and 

b. The choice of these parameters is important at high fre- 
quencies where unstable frequency ranges are found. At low 
frequency, covering radial orders up to n — 5, all modes are 
found to be stable for both sets of the a and b parameters. 
The values of | D p | are significantly higher than the values of 
~D P shown in the upper panel of Fig. 3. This clearly indi- 
cates that by neglecting the perturbed convective fluxes we 
ignore the dominant contribution to the damping. 

The results shown in Fig. 6 are applicable also to non- 
radial modes, because virtually all the contribution to D p 
arises in the upper layers where the value of £ has little in- 
fluence. However, damping effects in these layers have con- 
sequences in the deep interior. They change 7, and hence 
the amplitude behaviour in the g-mode propagation zone 
(see equation 3). We have seen in Section 3.2 (Figs. 2 and 
3) that ignoring driving effects in these layers reduces the 
trapping. Adding damping there would reduced it further. 
Larger inertiae imply lower amplitudes for stochastically ex- 
ited modes, and indeed we should not expect a detection of 
stochastically excited nonradial modes in giants. 
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Figure 6. Absolute values of energy dissipation rates, \D V \ (top) 
and energy generation rate PA 2 (bottom) (P is the rate of en- 
ergy injected into the modes by the fluctuating Reynolds stresses 
and A is defined by equation 22). The energy generation rates are 
expressed in units of the solar luminosity. Two sets for the con- 
vection parameters a and b were used. Thick curves in the upper 
panel indicated the frequency range where radial modes are found 
to be unstable (D p < 0). 



4.2 Amplitudes of stochastically driven radial 
modes 

The amplitudes of intrinsically stable stochastically driven 
radial modes were estimated in the manner of Houdek et al. 
(1999): 



Vs = 



(19) 



where here Pq is the noise generation rate injected into a 
mode through the fluctuating Reynolds stresses, the expres- 
sion for which we adopted from Balmforth (1992b) (see also 
Houdek etal. 1999). The damping rate is 77 = D p /2Iu> 2 — 
—7, and I w — IR~ 2 in our notation. For radial modes the 
total energy dissipation rate D is D p . The linear stability 
analysis also provides the parameter A, which is the ratio 
of the relative luminosity to the relative velocity amplitude, 
computed at the surface (i.e. outermost meshpoint) of the 
star. The bolometric relative luminosity amplitude then be- 
comes 



A 



SR 



A 



V a 



6L 

L " R "loR' 
and from equation (21) we obtain, 

5L - X 



(20) 



(21) 
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O Model (a, = & =300) 

▲ Model (a =900, f> =2000) 




Figure 7. Amplitudes and frequencies of oscillations in a UMa 
from Buzasi et al. (2000) compared with model calculations ob- 
tained for two sets of convection parameters. The width of the 
shadowed rectangles corresponds to the uncertainty of the fre- 
quency data. Calculated amplitudes for the n = 1, 2, and 3 modes 
are 0.5, 2.7, and 5.8 /imag, respectively. The observed peaks at 
v = 34.9 ± 0.6 and v = 43.6 ± 0.9 fiHz, which are above the 
calculated acoustic cut-off frequency, are not shown. 



where P = PqI u and I n = I/3MR 2 ; I n is the dimensionless 
modal inertia plotted in Fig. 3. In the lower panel of Fig. 6, 
we plot the quantity \ 2 P. All the quantities plotted in this 
figure are applicable also to nonradial modes of low degree. 
However, for nonradial modes we have to take into account 
the damping effects in the g-mode propagation zone. With 
the help of equation (21) and the data given in Fig. 3 we 
can evaluate amplitudes for radial modes with D p < 0. 

In Fig. 7, we compare radial-mode frequencies and am- 
plitudes calculated for M a with the observational data of 
a UMa. Keeping in mind the large observational errors and 
the fact that we have made no effort to adjust model pa- 
rameters to fit the frequencies, we regard the agreement of 
frequencies as satisfactory. On the other hand, the disagree- 
ment between the amplitudes is very serious: the observed 
amplitude at n = 1 exceeds the predicted value by three 
orders of magnitude, and the frequency dependence of the 
amplitudes differ drastically. 

An additional difficulty is presented by the presence of 
the two peaks above the acoustic cut-off frequency. Such 
high-frequency peaks are observed in the Sun, but with am- 
plitudes much lower than those below the acoustic cut-off. 
The two highest-frequency peaks in a UMa have amplitudes 
of about 0.2 mmag, which are similar to most of the other 
peaks. We should stress that the amplitude estimates in 
Fig. 7 were obtained using the pulsation modes of a model 
with an atmosphere based on model C of Vernazza, Avrett 
& Loeser (1981). That atmosphere has an acoustic cut-off 
frequency of 32.4/iHz at the temperature minimum, which is 
lower than the two highest frequencies of the observed peaks. 
A more realistic atmosphere might have a higher value than 
this, which itself is higher than the value for the Eddington 
grey atmosphere quoted in Table 2 and used by Guenther 
etal. (2000). 

The amplitudes of stochastically excited nonradial 
modes, including those that are most efficiently trapped in 
the acoustic cavity, are expected to have values much lower 
than those of corresponding radial modes. Equation (20) 
applies to nonradial modes if the contributions to I n and D 



from the g-mode propagation zone are included. The values 
of both I n and D g are substantially larger than those plotted 
in the lower panels of Fig. 2 and 3, owing to the damping 
effect in the outer layers. 



5 SOLAR-LIKE OR MIRA-LIKE 
OSCILLATIONS ? 

The results of the stability analysis presented in Section 4.1 
seem to exclude an interpretation of the low-frequency part 
of the a UMa oscillation spectrum in terms of self-excited 
modes. Indeed, the damping effect of convection exceeds by 
a large margin the driving effect of the opacity perturbation. 
However, there still seems to be a greater chance for an in- 
terpretation in terms of Mira-like excitation than in terms 
of solar-like excitation. The trend of calculated amplitudes 
is determined mainly by the factor I^ 1 ^ 2 in equation (21), 
which is the most reliably calculated quantity in the expres- 
sion. One may contemplate that for the first three modes D p 
is really much lower than what we calculated, but this op- 
tion would require near cancellation of damping and driving 
effects; it is more plausible that the quantity is less than zero 
and that the modes are unstable. The option that remains 
is an increase of X 2 E by four to six orders of magnitude. 

Whatever is the correct answer, the required changes 
are bound to be related to the way in which we treat the in- 
teraction between pulsation and convection. Our treatment, 
like most of those that have been used, is based on the 
mixing-length formalism, and we know that it is an inad- 
equate tool for describing the mean properties of convec- 
tion. In studies of acoustic mode damping and excitation we 
have to consider more detailed aspects of the dynamics of 
convection. The alternative is a hydrodynamical simulation. 
This has already been applied to solar radial oscillations by 
e.g. Stein and Nordlund (2001). It is to be hoped that be- 
fore long this approach will become applicable to red-giant 
oscillations too. 

How could future observational work on red-giant vari- 
ability help us? One possibility is the disproof of genuine 
pulsations in a UMa and in other red giants. Short-term 
variability could be a direct manifestation of convection, like 
large-scale granulation. This would result in progress being 
slow: compare how much we have learned in the past from 
the Sun's granulation with what we have learned from its 
oscillations. 

One very promising observational approach to the solar- 
like vs Mira-like alternatives is repeating the analysis of 
Cook's etal. (1997) MACHO data with a much longer time 
base or of extensive data from another microlensing projects 
such as OGLE (Udalski etal. 1997). Much improved fre- 
quency and amplitude resolution is expected. Showing that 
the ridges extend from a few days to hundreds of days with 
a continuous amplitude increase might strengthen the Mira- 
like interpretation. 

Other observational evidence supporting a Mira-like in- 
terpretation would be the identification of nonradial modes. 
We have seen in Section 3 that if radial modes are unsta- 
ble some i = 1 and 2 modes should be unstable too. If the 
modes are stable, then, as we discussed at the end of Sec- 
tion 4.2, nonradial modes will be excited stochastically, but 
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their amplitudes will be much lower than those of their ra- 
dial counterparts. 

It could be possible that low-order modes in a UMa are 
Mira-like whilst those of higher order are solar-like. This 
could also be the case for the multiperiodic M-type giants 
found by Koen and Laney (2000). In some of these stars 
the frequency ratio exceeds 10, and there is no doubt that 
the highest frequencies exceed the acoustic cut-off frequency. 
Two of the peaks in a UMa, as we have already noted, are 
also above the acoustic cut-off frequency, but that does not 
necessarily produce pulsational stability (cf. Balmforth et al. 
2001) 

Regardless of what the excitation mechanism is, the 
data on normal-mode frequencies will be very useful as a 
constraint on stellar parameters and models. Prospects of 
detecting nonradial modes is particularly interesting in this 
context. 
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